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Abstract 

In order to make quantitative predictions for jet cross sections in perturba- 
tive QCD, it is essential to calculate them to next-to-leading accuracy. This 
has traditionally been an extremely laborious process. Using a new formalism, 
imaginatively called the dipole formalism, we are able to construct a completely 
general algorithm for next-to-leading order calculations of arbitrary jet quanti- 
ties in arbitrary processes. In this paper we present the basic ideas behind the 
algorithm and illustrate them with a simple example. 
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1 Introduction 



Most of the recent progress in the understanding of strong interaction physics at large 
momentum transfer has been due to the comparison between precise experimental data 
and very accurate QCD calculations to higher perturbative orders |T|]. 

These higher-order computations have been carried out over a period of about fifteen 
years, often long after the accuracy of experimental data has made them necessary, be- 
cause of the difficulties in setting up a general and straightforward calculational procedure. 
The physical origin of these difficulties is in the necessity of factorizing the long- and 
short-distance components of the scattering processes and is reflected in the perturbative 
calculation by the presence of divergences. 

In general, when evaluating higher-order QCD cross sections, one has to consider real- 
emission contributions and virtual-loop corrections and one has to deal with different kind of 
singularities. The customary ultraviolet singularities, present in the virtual contributions, 
are removed by renormalization. The low-momentum (soft) and small-angle (coUinear) 
regions instead produce singularities both in the real and in the virtual contributions. In 
order to handle these divergences, the observable one is interested in has to be properly 
defined. It has to be a jet quantity, that is, a hadronic observable that turns out to be 
infrared safe and either coUinear safe or collinear factorizable: its actual value has to be 
independent of the number of soft and collinear particles in the final state (see Sect. ^ 
for a formal definition). In the case of jet quantities, the coherent sum over different 
(real and virtual) soft and collinear partonic configurations in the final state leads to the 
cancellation of soft singularities. The left-over collinear singularities are then factorized 
into the process independent structure and fragmentation functions of partons (parton 
distributions), leading to predictable scaling violations. As a result, jet cross sections 
are finite (calculable) at the partonic level order by order in perturbation theory. All the 
dependence on long-distance physics is either included in the parton distributions or in non- 
perturbative corrections that are suppressed by inverse powers of the (large) transferred 
momentum Q that controls the scattering process. 

Because of this pattern of singularities, QCD calculations of jet cross sections beyond 
leading order (LO) are very involved. Owing to the complicated phase space for multi- 
parton configurations, analytic calculations are in practice impossible for all but the sim- 
plest quantities, but the use of numerical methods is far from trivial because soft and 
collinear singularities present in the intermediate steps have first to be regularized. This 
is usually done by analytic continuation to a number of space-time dimensions d = A — 2e 
different from four, which greatly complicates the Lorentz algebra in the evaluation of the 
matrix elements and prevents a straightforward implementation of numerical integration 
techniques. Despite these difficulties, efficient computational techniques have been set up, 
at least to next-to- leading order (NLO), during the last few years. 

There are, broadly speaking, two types of algorithm used for NLO calculations: one 
based on the phase-space slicing method and the other based on the subtraction method^ 
The main difference between these algorithms and the standard procedures of analytic cal- 
culations is that only a minimal part of the full calculation is treated analytically, namely 
only those contributions giving rise to the singularities. Moreover, for any given pro- 

'''We refer the reader to the Introduction of Ref. Q for an elementary description of the basic difference 
between the two methods. 
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cess, these contributions are computed in a manner that is independent of the particular 
jet observable considered. Once every singular term has been isolated and the cancella- 
tion/factorization of divergences achieved, one can perform the remaining part of the cal- 
culation in four space-time dimensions. Although, when possible, one still has the freedom 
of completing the calculation analytically, at this point the use of numerical integration 
techniques (typically, Monte Carlo methods) is certainly more convenient. First of all, 
the numerical approach allows one to calculate any number and any type of observable 
simultaneously by simply histogramming the appropriate quantities, rather than having to 
make a separate analytic calculation for each observable. Furthermore, using the numerical 
approach, it is easy to implement different experimental conditions, for example, detector 
acceptances and experimental cuts. In other words, the phase-space slicing and subtrac- 
tion algorithms provide the basis for setting up a general-purpose Monte Carlo program 
for carrying out arbitrary NLO QCD calculations in a given process. 

Both the slicing and the subtraction ^ methods were first used in the context of 
NLO calculations of three-jet cross sections in e~^e~ annihilation. Then they have been 
applied to other cross sections adapting the method each time to the particular process. 
Only recently has it become clear that both algorithms are generalizable in a process- 
independent manner. The key observation is that the singular parts of the QCD matrix 
elements for real emission can be singled out in a general way by using the factorization 
properties of soft and coUinear radiation [Q. 

At present, a general slicing algorithm is available for calculating NLO cross sections 
for any number of jets both in lepton ||^ and hadron [j^ collisions. To our knowledge, 
fragmentation processes have been considered only in the particular case of direct-photon 
production [§]. The complete generalization of this method to include fragmentation func- 
tions and heavy flavours is still missing. 

As for the subtraction algorithm, a general NLO formalism has been set up for comput- 
ing three-jet observables in e~^e~ annihilation and cross sections up to two final-state 
jets in hadron collisions^. Also the treatment of massive partons has been considered 



in the particular case of heavy-quark correlations in hadron collisions [12 



In this paper, we present the basic idea to set up a completely general version of the 
subtraction algorithm. This generality is obtained by fully exploiting the factorization 
properties of soft and coUinear emission and, thus, deriving new improved factorization 
formulae, called dipole factorization formulae. They allow us to introduce a set of universal 
counter-terms that can be used for any NLO QCD calculation. 

For the purpose of illustration, in this short contribution we describe the implementation 
of our general method to the calculation of jet cross sections in processes with no initial- 
state hadrons, typically e~^e~ annihilation. Full details of the method and its application 



to all the other hard-scattering processes will appear elsewhere |13 



We first recall the main features of the subtraction method in Section |^. Then in 
Section |^ we present our dipole factorization formulae, for the case in which there are no 
incoming QCD partons. These allow us to calculate cross-sections for an arbitrary number 
of jets in e~^e~ annihilation, which we do in Section ^ After briefly recapping the resulting 
formulae in Section ^ we illustrate them with a simple example in Section ^, e'^e' — > 3 
jets. Finally in Section ^ we give a summary and outlook. 

^An extension of the method for three-jet cross sections has been recently presented [Q. 
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2 The subtraction procedure 



Suppose we want to compute a jet cross section a to NLO, namely 

^LO I „NLO (-\\ 

Here the LO cross section is obtained by integrating the exclusive cross section da in 
the Born approximation over the phase space for the corresponding jet quantity. Suppose 
also that this LO calculation involves m partons in the final state. Thus, we write 

a^o = f da"" , (2) 



where, in general, all the quantities (QCD matrix elements and phase space) are evaluated 
in (i = 4 — 2e space-time dimensions. However, by definition, at this LO the phase space 
integration in Eq. (§) is finite so that the whole calculation can be carried out (analytically 
or numerically) in four dimensions. 

Now we go to NLO. We have to consider the exclusive cross section da^ with m + 1 
partons in the final-state and the one-loop correction da^ to the process with m partons 
in the final state: 

J Jm+1 Jm 

The two integrals on the right-hand side of Eq. (Q) are separately divergent if c? = 4, 
although their sum is finite. Therefore, before any numerical calculation can be attempted, 
the separate pieces have to be regularized. Using dimensional regularization, the diver- 
gences (arising out of the integration) are replaced by double (soft and collinear) poles 
and single (soft, collinear or ultraviolet) poles 1/e. Suppose that one has already carried out 
the renormalization procedure in da^ so that all its ultraviolet poles have been removed. 

The general idea of the subtraction method for writing a general-purpose Monte Carlo 
program is to use the identity 



da 



NLO 



da^ - da' 



+ da'' + da'' , (4) 



where da^ is a proper approximation of da^ such as to have the same pointwise singular 
behaviour (in d dimensions) as da^ itself. Thus, da^ acts as a local counterterm for da^ 
and, introducing the phase space integration. 



^NLO 



m+1 



da^ - da' 



+ da''+ da"" , (5) 

J m+l Jm 



one can safely perform the limit e ^ under the integral sign in the first term on the 
right-hand side of Eq. (^. Hence, this first term can be integrated numerically in four 
dimensions. 

All the singularities are now associated to the last two terms on the right-hand side of 
Eq. (^. If one is able to carry out analytically the integration of da"^ over the one-parton 
subspace leading to the e poles, one can combine these poles with those in da^ , thus 
cancelling all the divergences, performing the limit e ^ and carrying out numerically the 
remaining integration over the m-parton phase space. The final structure of the calculation 
is as follows 

(6) 



^NLO 



[ Ida^^o - daf^ol + f da^ + f da 
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and can be easily implemented in a 'partonic Monte Carlo' program, which generates 
appropriately weighted partonic events with m + 1 final-state partons and events with m 
partons. 

The key for the subtraction procedure to work is obviously the actual form of da^. One 
needs to find an expression for da^ that fulfils the following properties: i) for any given 
process, da"^ has to be obtained in a way that is independent of the particular jet observable 
considered; ii) it has to exactly match the singular behaviour of da^ in d dimensions; Hi) its 
form has to be particularly convenient for Monte Carlo integration techniques; iv) it has to 
be exactly integrable analytically in d dimensions over the single-parton subspaces leading 
to soft and collinear divergences. 

In Ref. [0, a suitable expression for da^ for the process e+e~— >• 3 jets was obtained by 
starting from the explicit expression (in d dimensions) of the corresponding da^ and by 
performing extensive partial fractioning of the 3 + 1-parton matrix elements, so that each 
divergent piece could be extracted. This is an extremely laborious and ungeneralizable 
task, in the sense that having done it for e^e^^ 3 jets does not help us to do this for, say, 
e"'"e~— > 4 jets or for any other process. 

In Ref. [0], the general properties of soft and collinear emission were first used (in the 
context of the subtraction method) to construct da^, for one- and two-jet production in 
hadron collisions, in a way that is independent of the detailed form of the corresponding 

The central proposal of our version of the subtraction method is that one can give a 
recipe for constructing da"^ that is completely process independent (and not simply indepen- 
dent of the jet observable). Starting from our physical knowledge of how the m + 1-parton 
matrix elements behave in the soft and collinear limits that produce the divergences, we 
introduce improved factorization formulae, called dipole formulae (see Sect. §), which allow 
us to obtain in a straightforward way (see Sect. ^ a counter-term da^ satisfying all the 
properties listed above. 

3 Dipole factorization formulae 

Notation 

In general we use dimensional regularization in ci = 4 — 2e space-time dimensions and 
consider d — 2 helicity states for gluons and 2 helicity states for massless quarks. This 
defines the usual dimensional-regularization scheme. Other dimensional-regularization pre- 
scriptions can be used, as well [ |T3|1 . 

The dimensional-regularization scale, which appears in the calculation of the matrix 
elements, is denoted by /i. In the perturbative calculation of physical cross sections, after 
having combined the renormalized matrix elements, the dependence on fi exactly cancels 
and is replaced by the dependence on the renormalization scale ^r. Therefore, in order to 
avoid a cumbersome notation, we set /i = fiR. 

Throughout the paper, as stands for Q;s(/i), the NLO QCD running coupling evaluated 
at the renormalization scale /z. The actual value of the QCD coupling as{fi) depends on 
the renormalization scheme used to subtract the ultraviolet divergences from the (bare) 
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one-loop matrix element (or, equivalently, from da^ in Eq. (|^)). 

The d-dimensional phase space, which involves the integration over the momenta {pi, 
Pro] of m final-state partons, will be denoted as follows 



{2nY6^''\pi + ...+p^-Q) = d(Pm{pi,.-,Pm;Q) . (7) 



In the case of processes without initial-state QCD partons (e+e~-type processes), the 
{tree-level) matrix element with m QCD partons in the final state has the following general 
structure (non-QCD partons, namely 7*, Z^, W"^, . . ., carrying a total incoming momentum 
Qfj,, are always understood) 

M'^-'--^'^-'-ip,,...,p^) (8) 

where {ci,...,Cm}, {si,...,Sm} and {pi,..,pm} are respectively colour indices (a = 1,..., 
A^^ — 1 different colours for each gluon, a = 1, ..,Nc different colours for each quark or an- 
tiquark), spin indices {fi = 1, ...,d for gluons, s = 1, 2 for massless fermions) and momenta. 

It is useful to introduce a basis {|ci, Cm > ®|si, >} in colour + hehcity space 
in such a way that 

A<^'-'''"''''-''"(Pl,-.,Pm) = (< Ci,...,Cm\® < Si,...,S^|)|l,...,m >^ . (9) 

Thus |l,...,m >jn is a vector in colour + helicity space. According to this notation, the 
matrix element squared (summed over final-state colours and spins) |A^„p can be written 
as 

l^mT = m< 1, •••,"^|l,.-,"^>m • (10) 

As for the colour structure^, it is convenient to associate a colour charge Tj with the 
emission of a gluon from each parton i. We thus define the square of colour- correlated 
tree-amplitudes as follows 

\Mli\^ = l,...,m|T, ■Tfc|l,...,m>^ 

where T^^, = ifcat (colour-charge matrix in the adjoint representation) if the emitting par- 
ticle i is a gluon and T^^ = (colour-charge matrix in the fundamental representation) if 
the emitting particle z is a quark (in the case of an emitting antiquark T^^ = = —f^^). 
It is straightforward to check that the colour-charge algebra is: 

T,-Tj = Tj-T, ifz^j; = C,, (12) 

where Cj is the Casimir operator, that is, Ci = Ca = A^'c if i is a gluon and Ci = Cp = 
{N^ — l)/2Nc if i is a quark or antiquark. 

In this notation, each vector |1, ...,m >m is a colour singlet, so colour conservation is 
simply 

m 

5]T,|l,...,m>„=0. (13) 

i=l 



^Within our formalism, there is no need to consider the decomposition of the matrix elements into 
colour subamplitudes, as in ||^,p^. 
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Dipole formulae 

The real contribution da^ to the NLO cross section in Eq. @ is proportional to the tree- 
level matrix element A4m+i for producing m + 1 partons in the final state. The dependence 
of |7Wm+iP on the momentum pj of a final-state parton j is singular in two different phase- 
space regions: when the momentum pj vanishes {soft region) and/or when it becomes 
parallel to the momentum pi of another parton in Aim+i {coUinear region). This singular 
behaviour of is well-known Pjr5| and universal. Indeed, in the soft and collinear 

limits, M.m+1 is essentially factorizable with respect to M-m-, the tree-level amplitude with 
m partons, and the singular factor only depends on the momenta and quantum numbers 
of the QCD partons in A^^- 

We thus introduce the following dipole factorization formula: 

\Mm+l{pi, ...,Pm+l)? =m+l< 1, 1||1, ...,m+ 1 >m+l= '^ij,k{Pl^ ■■■,Pm+l) + ■■■ 

(14) 

where . . . stands for terms that are not singular in the limit Pi ■ pj ^ (i.e. when i and j 
become collinear or when either i or j is soft) and the dipole contribution 'Dij^k is given by 

1^ij,k {Pl, ■■■,Pm+l) = ^ — (15) 

2pi ■ Pj 

~ ~ Tj^ ■ T - ■ ~ ~ 

■ m< 1) ■■I'^j 1 ..yTn + 1| —2 ij.k 1 1, . . , , . . , fc, . . , m -|- 1 >m . 

The m-parton matrix element on the right-hand side of Eq. (|T3|) is obtained from the 
original m + 1-parton matrix element by replacing a) the partons i and j with a single 
parton ij (which plays the role of emitter) and h) the parton k with the parton k (which 
plays the role of spectator). All of the quantum numbers (colour, flavour) except momenta 
are assigned as follows. The spectator parton k has the same quantum numbers as k. The 
quantum numbers of the emitter parton ij are obtained according to their conservation 
in the collinear splitting process ij —>■ i + j (i.e. 'anything + gluon' gives 'anything' and 
'quark + antiquark' gives 'gluon'). 

The momenta of the emitter and the spectator are deflned as follows 

Pij =Pi+ Pj - z-^^^Pk > Pk = -r^- Pk ' (16) 

-L yij,k ^ Vijfk 

where the dimensionless variable Uij^k is given by 

PiPj /^7^ 
yij,k = ^ ^ • 17 

PiPj + PjPk + PkPi 

In the bra-ket on the right-hand side of Eq. (|l^), Tij and Tk are the colour charges of 
the emitter and the spectator and Vij^k are matrices in the helicity space of the emitter. 
These matrices, which depend on i/ij^k and on the kinematic variables Zi, zj: 

~ _ PiPk _ PiPk - _ PjPk _ PjPk _~ /.ON 



PjPk + PiPk PijPk PjPk + PiPk PijPk 
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are universal factors related to the li-dimensional Altarelli-Parisi splitting functions [15]. 
For fermion + gluon splitting we have (s and s' are the spin indices of the fermion ij in 
< and > respectively) 



< s\Vq^g^^k{zi;yijM)\s' > = Snfi ^as Cp 



1- Zi{l- yij^k) 



(1 + Zi) - e(l - Zi) 



^qigj,k ^ss' ■ 



(19) 



For quark + antiquark and gluon + gluon splitting we have (/i and u are the spin indices 
of the gluon ij in < ..\ and .. > respectively) 



< ^^\^q^qJ,k{^i)W >= 8^/^ '"S Tr 



PiPj 



(20) 



< IJ'\Vg,g^,k{zi;yij,k)W >= Ifivr/x 'as Ca 
+ 



1- Zi{l- yij^k) 



7 - 2 ) + (1 - e)^ {~z,pf - ~z,p^) {~z,pt - z,p'<) 

^ - Zj[l - yij,k) ) PiPj ^ 



y^.M- (21) 



The factorization formula in Eq. (|T^) has a dipole structure with respect to the colour 
and spin indices of the factorized partons. As shown in Ref. |r3|, Eqs. (n.T^) coincide with 
the soft-gluon ||^ and Altarelli-Parisi []15[ factorization formulae respectively in the soft 
and coUinear limits. However, Eqs. (0,|lD are completely well-defined also outside these 



limiting regions of the phase space. Indeed, in the factorized m-parton matrix element 

0) and, performing the 



pi 



both the emitter ij and the spectator k are on-shell {pf 
replacement A;} — {zj,A;}, momentum conservation is implemented exactly: 



Pi +P^j + Pk 



Ptj + Pk 



(22) 



The importance of these kinematical features is twofold. Firstly, momentum conservation 
leads to a smooth interpolation between the soft and collinear limits, thus avoiding double 
counting of overlapping soft and collinear singularities. Secondly, the definition (|TBp of 
the dipole momenta allows us to factorize exactly the m + 1-parton phase space into an 
m-parton subspace times a single-parton contribution (see Eqs. ( PT|J5^ )). The first property 
allows us to construct a counter-term da"^ that produces a pointwise cancellation of the 
singularities of da^ as in Eq. (|^). The second property makes this counter-term fully 
integrable analytically over the subspace leading to soft and collinear divergences. 



4 The calculation of jet cross sections 

The dipole formulae form the basis for our general algorithm for NLO jet calculations, as we 
describe below. First we define the leading order cross section and give a formal definition 
of the requirements a jet definition must fulfil. Then we introduce the subtraction term, 
which cancels all singularities of the real matrix element, and show how it can be integrated 
in d dimensions to cancel the singularities of the virtual matrix element. 
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Leading order and jet definition 

The Born-level cross section in Eq. (|]) has the following expression 

da^ = d(j)rn{pi, ■■■,Prn]Q) \Mrn{pl, ■■■,Pm)\^ F^"^\pi, ■■■,Pm) , (23) 

{m} "^{"i} 

where N'm includes all the factors that are QCD independent, J2{m.} denotes the sum over 
all the configurations with m partons, dcpm is the partonic phase space in Eq. (|^), S^m} is 
the Bose symmetry factor for identical partons in the final state and Aim is the tree-level 
matrix element. 

The function Fj™"-* (^p^^ ...,Pm) defines the jet observable in terms of the momenta of the m 
final-state partons. In general, Fj may contain ^-functions (thus, Eq. ( P^D defines precisely 
a cross section), 5-functions (Eq. (p3D defines a differential cross section), numerical and 
kinematical factors (Eq. ( p3D refers to an inclusive observable), or any combination of 
these. The essential property of Fj™-* is that the jet observable we are interested in has to 
be infrared and collinear safe. That is, it has to be experimentally (theoretically) defined in 
such a way that its actual value is independent of the number of soft and collinear hadrons 
(partons) produced in the final state. In particular, this value has to be the same in a 
given m-parton configuration and in all m + 1-parton configurations that are kinematically 
degenerate with it (i.e. which are obtained from the m-parton configuration by adding a 
soft parton or replacing a parton with a pair of collinear partons carrying the same total 
momentum). These properties can be simply restated in a formal way, as follows 

^i"^'HPl,-,P.,-,Pn+l)^i^i"^(Pl,...,Pn+l) if P,-^0 , (24) 

Fj'^^\pi,-;Pt,-,Pj,--,Pn+i) ^ F!f\pi,..,pi+pj,..,pn+i) if PiWPj, (25) 
)(pi, ...,p„) ^0 if ■ p,-^ . (26) 

Equations (0) and (^) respectively guarantee that the jet observable is infrared and 
collinear safe for any number n of final-state partons, i.e. to any order in QCD perturbation 
theory. Equation (^) defines the LO cross section, that is, it ensures that the Born-level 
cross section da^ in Eq. (^3]) is well-defined (i.e. finite after integration) in c? = 4 dimensions. 

Next-to-leading order: the subtraction term 

The real contribution da^ to the NLO cross section in Eq. (|^) has the same expression 
as da^ in Eq. (^), apart from the replacement m m + 1. In particular, the m-parton 
matrix element Aim is replaced by Aim+i- Therefore an explicit and general form for the 
local counter-term da^ in Eq. (^ is provided by the dipole factorization formula (pjlf): 

(icr^ = Min d(f)m+l{pi, ■■■,Pm+i;Q) 

{m+1} 'J{m+1} 

• Y '^ijAPl^---^Pni+l) Flj"'\pi,..pij,Pk,..,Pm+l) ■ (27) 

pairs k^i,j 

Here T>ij k{pi, ...,pm+i) is the dipole contribution in Eq. (P^) and Fj^\pi, ..pij,pk, ..,pm+i) 
is the jet function for the corresponding m-parton state {pi, ..pij,pk, ..,Pm+i 

}. Note that 
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this is completely independent of Pi, which is how we are able to integrate da^ analytically 
over the phase-space of i without any information about the form of Fj, as we perform 
below. 

We can check that the definition (^) makes the difference {dcr^ — da"^) integrable in 
d = 4 dimensions. Its explicit expression is 

da^-da"^ = Mn Y d4>m+i{pi, ■■■,Pm+i; Q) — 

{m+l} 



\Mm+l{pi, F'f'^^\pi, 
Y 1^tj,kiPl^--^Pm+l) Fi"'\pi,--Pij,Pk,-,Pm+l)\ ■ (2^ 



pairs k^i,j 
iJ 

Each term in the curly bracket is separately singular in the soft and collinear regions. 
However, as stated in Sect. in each of these regions both the matrix element M.m+i and 
the phase space for the m + 1-parton configuration behave as the corresponding dipole 
contribution and dipole phase space: 

\Mm+l{Pl,--,Pm+l)\^ Vij^kiPl, ■■■,Pm+l) , (29) 

{pi,..Pi,..Pj,..Pk,..,Pm+l} {Pl,-Pij,Pk,-;Pm+l} ■ (30) 



Thus, because of Eqs. (^) and (p5D, the singularities of the first term in the curly bracket 
are cancelled by similar singularities due to the second term. On the other hand, each 
dipole 'Dij k in Eq. ( [T3| ) has no other singularities but those due to the m-parton matrix 
element \1, , .., k, ..,m + 1 >„. Because of Eq. (^6]) , these singularities are screened 
(regularized) by the jet function Fj^\pi, ..pij,pk, ..,Pm+i) in the curly bracket of Eq. (pSf ). 

Note that this cancellation mechanism is completely independent of the actual form of 
the jet defining function but it is essential that da^ and da"^ are proportional to Fj"^^^^ 
and Fj"^^ respectively. Nonetheless, because the terms on the left- and right-hand sides of 
Eq.(p!^ depend on the same kinematic variables, {pi, ...,pm+i}, both da^ and da^ live on 
the same m + 1-parton phase space. Thus the numerical integration (in c? = 4 dimensions) 
of Eq. ( PSP via Monte Carlo techniques is straightforward. One simply generates an m + 1- 
parton configuration and gives it a positive (+ |7Vlm+iP) or negative {—Vij^k) weight. The 
role of the two different jet functions i^j™^^-* and Fj™^ is that of binning these weighted 
events into different bins of the jet observable. Any time that the generated configuration 
approaches a singular region, these two bins coincide and the cancellation of the large 
positive and negative weights takes place. 

Note, also, that the helicity dependence of the splitting kernels Vij^k in Eqs. (|^,^) is 
essential if da^ is to act as a local counterterm that makes [da^ — da^] integrable in four 
dimensions. Indeed, the parton azimuthal correlations due to this dependence are not only 
essential in the most general case when Fj explicitly depends on them, but even when it 
does not]]]. 

^In this case the evaluation of J^^^ da^ in four dimensions usuaUy involves double angular integrals of 

the type J^^ dcos9 J^^ dip cosiy9/(l — cos6'), where (p is the azimuthal angle. These integrals are mathe- 
matically ill-defined. If their numerical integration is attempted, one can obtain any answer whatsoever, 
depending on the detail of the integration procedure. Performing the integral analytically before going to 
4 dimensions, one obtains J^^ dcos9 J^^ dip cos(/?/(l — cos 6*) sin^^*^ ^sin"^"^ tp = 0. 
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Next-to-leading order: integral of the subtraction term 

Having discussed the four- dimensional integrability of {da^ — da^), the only other step 
we have to consider is the dimensional analytical integrability of da^ over the one-parton 
subspace leading to soft and collinear divergences. In this respect, our dipole formalism is 
particularly efficient and simple. The definition ([T6|) of the dipole momenta allows us to 
exactly factorize the phase space of the partons i,j, k into the dipole phase space times a 
single-parton contribution, as follows 



{pi,..,Pi,Pj,Pk,.., Pm+1 ]Q) = d(f)m {pi,..,Pij,Pk,.., Pm+1 ] Q) [dpi {pij , Pfc)] 



(31) 



where 



[dpi{.Pij,Pk)] 



d'^Pi 



^'5+(r')0(i-5.)0(i-2/..,a 



(32) 



and the kinematic variables yij^k and Zi are defined in Eqs. (p!7| , p!8D . 

Inserting Eq. (|3T|) and the explicit expression (|15]) for Vij^k into Eq. (^), we can factorize 
completely the Pi dependence and carry out the integration over the phase space region (^) 
with the m parton momenta {pi, ..,pij,pk, ..,Pm+i} kept fixed. Remarkably, this integration 
can be exactly performed in closed analytical form in any number of space-time dimensions. 
As shown in detail in Ref. [|1^], after integration the spin correlations between Vij^k and 
|1, .., ij, .., fc, .., m + 1 >m vanish and only colour correlations survive. The final result for 
Jm+i dc^ can be written in terms of an m-parton integral of the LO (colour-correlated) 



matrix element times a factor 13 



/ da^= f 






f 


'm+l Jm 






Jm 



/ N'inY. d4>m{Pl,---,Pm]Q) 

J m r 1 



1 



{m} "^i™} 

• ^< 1, ...,m| J(e) |1, ...,m >^ Fj"'^(pi, ...,p^) . (33) 

Comparing Eqs. (^) and (^3]), we see that the integration of da^ over the one-parton 
subspace that produces soft and collinear singularities leads to an expression completely 
analogous to da^ . One should simply replace the matrix element squared = 



i< 1, m|l, m >m in da with 



1 m\ lie) 11 m >, 



(34) 



where 7(e) is an insertion operator that depends on the colour charges and momenta of the 
m final-state partons. Its explicit expression is ||13||: 



J(pi,...,p^;e) 



as 



where the singular factors Vj(e) have the following e-expansior||| 



(35) 



VAt) 



vr 



+ -1i- + -1^ + Ki + C(e) 



(36) 



"Their exact expressions in any number d = 4 — 2e of dimensions are given in |T^]. 
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with {Tfi = 1/2 and Nj is the number of flavours) 

%=g,?=|Ci., ^^=g = fCA-lTnNf, (37) 

In the calculation of the NLO cross section Eq. (|33|) has to be combined with 
the virtual contribution, whose expression in terms of the (renormalized) one-loop matrix 
element is the following 

da^ =J\finY^ d(t)m{Pl,-;Pm]Q) I A^m(Pl , • • • , Pm) | (l_ioop) Fj^\pi^ ■■■,Pm) ■ (39) 

{m} '-'{'»} 

As discussed in Ref. |jl3|, the addition of these two contributions correctly produces the 



cancellation of all the e-poles, thus leading to a flnite NLO cross section. 



5 Final formulae 

The flnal results of the application of our algorithm to the calculation of jet cross sections 
with no hadron in the initial state are summarized below. 

The full QCD cross section in Eq. (|l|) contains a LO and a NLO component. Assuming 
that the LO calculation involves m flnal-state partons, the LO cross section is given by 



/ rfa^= /d<l>(-) |M^(pi,...,pJpFf)(pi,...,pJ , (40) 



where M.m is the tree-level QCD matrix element for producing m partons in the flnal 
state and the function Fj"*^ deflnes the particular jet observable we are interested in (see 
Eqs. ( P^P^D for the general properties that Fj^^ has to fulfll). The factor ci^^*") collects 



all the relevant phase space factors, i.e. all the remaining terms on the right-hand side 
of Eq. (P^. The whole calculation (phase space integration and evaluation of the matrix 
element) can be carried out in four space-time dimensions. 

According to the subtraction formula, Eq. (^), the NLO cross section is split into two 
terms with m + 1-parton and m-parton kinematics, respectively. The contribution with 
m + 1-parton kinematics is the following 

/ [dat,-dato] = /rf^^-^+^^l |-M^+i(pi,...,P™+i)rFf~''^(pi,...,p™+i) (41) 

~11 J2 '^ijAPl^-'^Pm+l) Fi"'\pi,--Pij,Pk,-;Pm+l) \ , 

pairs kj^ij 

where the term in the curly bracket is exactly the same as that in Eq. (pS]): Aim+i is the 
tree-level matrix element, Vij k is the dipole factor in Eq. (|15D and Fj"^^ is the jet deflning 
function for the corresponding m-parton state (note, again, the difference between the two 
jet functions i^j™^^-* and Fj^^ in the curly bracket). In spite of their original ci-dimensional 
deflnition, at this stage the full calculation is carried out in four dimensions. 
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The NLO contribution with m-parton kinematics is given by 



/ \da^ + f c/a^l (42) 

= J d^^""^ {\MmiPU ...,Pm)\l^loap) + m< I, -^ml I{e) 1 1, . . . , m }^^^ F]""^ (pi , . . . , p^) . 



The first term in the curly bracket is the one-loop renormalized matrix element squared 
for producing m final-state partons. The second term is obtained by inserting the colour- 
charge operator of Eq. (^5]) into the tree-level matrix element for producing m partons as 
in Eq. (0). These two terms have to be first evaluated in d = 4 — 2e dimensions. Then 
one has to carry out their expansion in e-poles (the expansion for the singular factors Vj(e) 
is given in Eq. (0)), cancel analytically (by trivial addition) the poles and perform the 
limit e — > 0. At this point the phase-space integration is carried out in four space-time 
dimensions. 



6 e+e — ^ 3 jets 

In this Section we consider the simplest non-trivial application of our algorithm, namely 
the calculation of three-jet observables in e'^e~ annihilation. Thus our formalism can be 
directly compared with that in Ref. Q]. 

The LO partonic process to be considered is e~^e~ qipi) + ^(^2) + oiPs)- The corre- 
sponding tree-level matrix element is denoted by A^3(pi, P2; Ps)- We use customary notation 
for the kinematic variables: is the square of the centre-of-mass energy, yy = 2pi -pj/Q"^ 
and Xi = 2pi ■ Q/Q"^. 

At NLO, two different real-emission subprocesses contribute: a) e~^e~ — > q{pi) + q{p2) + 
Qipz) + gipi)] b) e^e~ — >• q{pi) + q{p2) + q{p^) + ^(^4)- In addition one has to compute the 
one-loop correction to the LO process. 



The calculation of the subtracted cross section (41) for the subprocess a) involves the 
evaluation of the following dipole contributions: 2^13,2, ^^13,4, ^^14,2, ^^14,3, "^23,1, ^^23,4, "^24,1, 
^^24,35^^34,15^^34,2- The associatcd colour algebra can be easily performed in closed form 
because the several colour projections of the three-parton matrix element fully factorize[2[ 
Thus we do not need to compute any colour-correlated tree amplitudes and we find 

1 / C \ 

^^13,2(Pl5P25P35P4) = (1 - ^9193,2 l-M3(Pl35P25P4)r 5 



2piP3 V 2C, 
1 Ca 

2piP3 2Cp 

1 1 

2p3P4 2 



1^13a{Pi,P2,P3,P4) = TTTT" 7T7^ ^9193,4 1^3(^135 P25P4)P 5 (43) 
1^34,l{Pl,P2,P3,Pi) = T;:;^ ^ ^3334,1 ^^.(Pl 5 P2 5^34 ) ■ 



The dipole contributions "^23,15 ^^23,45 ^^34,2 are obtained respectively from ^^13,25 ^^13,45 "^34,1 
by the replacement pi ^ P2- Analogously, one can obtain 1^14,2 and 1^14,3 respectively from 
1)13 2 and 2^13,4 by the replacement p^ ^ p^, and ©24,1 and 1^24,3 respectively from 1^13,2 
and 1^13,4 by the replacement pi ^ P2, P3 ^ P4,- 

**If the LO matrix element involves two or three partons, the colour algebra can always be carried out 
in closed (factorized) form. 
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In the case of the subprocess b) we have to consider the following dipole contributions: 
^^12,3, ^^12,4, ^^14,2, "^^14,3, ^^23,1, ^^23,4, ^^34,1, ^^34,2- Performing the colour algebra we get 



^^34,l(pi,P2,P3,P4) 



11 

o K354,l^M^(Pl'P2,P34) 



2p3P4 2 



(44) 



and all the other dipoles are obtained by the corresponding permutation of the parton 
momenta. 

The splitting functions Vij^k of Eqs. (143| , |4^) are given in Eqs. ([19| - |2TD . The tensor 
T^u is the squared amplitude for the LO process e"'"e~ — >• qqg not summed over the gluon 
polarizations {fi and u are the gluon spin indices and —g^^T^^, = |A^3p). This can be 
easily calculated. In the case of jet observables averaged over the directions of the incoming 
leptons (un-oriented events) we find (in c? = 4 dimensions) 



^^^^(^1,^2,^3) 



1 "I" 



(45) 



where 



rp/lU 



r, PlP2 ^ _ ^ 1 - Xi p1p'( _ 2 1 - X2 P2P2 



I-X2 1-xi 



1 — Xi — X2 + X' 



X2 



+ (1 + + 1 



P1P3 pm_ 



X2 — X\ ~\- X-]^ 
1 — Xi 



1^ u fJ- u 

P2P3 P3P2 

+ Q2 



X2 — xi — X2j g 



(46) 



To complete the NLO calculation we also need the virtual cross section. In the case of 
un-oriented events, we take the one-loop matrix element in the MS renormalization scheme 
from Ref. (we use slightly different notation): 



\M3ipi,P2,P3)\ 



2 I A ^ / M2 1 / 47ryu^ 

(l-;oop) = l-^3(Pl,P2,P3)| ^ r(l-6) \ 

11 



1/23 



--[3Cf + — C 



6 



TT 



+ —{2Cf + Ca) 



as 
2tt 



[F{yi2,yi3,y23) + 0{e)] 



(47) 



where -^(1/12, ^/is, I/23) is defined in Eq. (2.21) of Ref. 



The explicit evaluation of the insertion operator /(e) in Eqs. (|34| , p5D gives: 



3< 1,2,3| /(e) |1,2,3 >3= \M3{pi,p2,P3)\ 



2 "s 



27rrfl 



47r/i^ 



{2Cf - CA)yii + Ca [y^i + fe'jj + - (27, + 7^ 



{2Cf - Ca) lnyi2 + -C^ ln(?/i3?/23) 



-7gln(?/i3y23) 



TT 



-y (2Ci. + Ca) + 2(7g + K,) +^g + Kg + 0{e) 



(4J 



13 



THRUST 



C-PARAMETER 






Figure 1: Coefficient of (as/27r)^ for the thrust and C-parameter distributions. The dotted 
histograms show the size of the statistical errors. 



Combining the one-loop matrix element (^7p with the result (^) according to Eq. and 
using the explicit expressions (|37|j3^) for 7^ and Ki, all the pole terms cancel. Note that as 
well as the pole terms, the closely related vr^ and In^ terms cancel: 



—loop) 

+ 3< 1,2,3| J(e) |1,2,3> 
3,_„ „ . . 1 



|^3(Pl,P2,P3)P 



27r 



2 {2Cf - Ca) In yi2 - -(5Ca - TnNf) ln(yi32/23) 



+2Cf + ^Ca - ^T^Nf 
9 9 



+ ^[F(|/i2,l/l3,y23)+C?(e)]] • 



(49) 



We have implemented these results as a working Monte Carlo programf^and the results 
are in good agreement with Ref. for all distributions shown there. As an example we 
show the NLO coefficients for the thrust and C-parameter distributions in Fig. |I|. We find 
that in general, the numerical convergence is similar to the program of Ref. |p, except 
close to the two-jet region in which ours becomes progressively better. More details of the 
Monte Carlo program and the generalization to oriented three-jet events will be presented 
elsewhere. 



7 Summary and outlook 

In this letter we have presented the basic idea to set up a completely general algorithm for 
calculating jet cross sections in NLO QCD. By general we mean that the algorithm applies 
to any jet observable in a given scattering process as well as to any hard-scattering process. 
The algorithm overcomes all the analytical difficulties related to the treatment of soft and 
collinear divergences in the perturbative expansion. The output of the algorithm is given 
in terms of effective matrix elements (the contributions in the curly bracket of Eqs. ( ^ , ^21) ) 
with built-in cancellation of soft and collinear singularities. These effective matrix elements 
can be numerically or analytically (whenever possible) integrated over the available phase 
space to compute the actual value of the NLO cross section. If the numerical approach is 

ttThe program can be obtained from http : //suryall . cern. ch/users/seymour/nlo/. 
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chosen, Monte Carlo integration techniques can be easily implemented to provide a general- 
purpose Monte Carlo program for carrying out NLO QCD calculations in any given process. 

Starting from the very general idea of the subtraction method, we have discussed how 
one can automatically construct a pointwise and integrable counter-term da^ for the real 
contribution da^ to the NLO cross section. This counter-term is independent of the details 
of the process under study and can be handled once and for all to define the effective 
matrix elements mentioned above. We have shown how this procedure works in practice 
for processes with no hadron in the initial state (typically, e~^e~ annihilation) and any 
number of jets in the final state. In this case, the final result is given by the master 



formulae in Eqs. (^,|^). We have provided explicit expressions for both the universal 
dipole factors Vij^k and /. Having these factors at our disposal, the only other ingredients 
necessary for the full NLO calculation, are the following 

• a set of independent colour projection^ of the matrix element squared at the Born 
level, summed over parton polarizations, in d dimensions; 

• the one-loop matrix element in d dimensions; 

• an additional projection of the Born level matrix element over the helicity of each 
external gluon in four dimensions; 

• the tree-level NLO matrix element in four dimensions. 



These few ingredients are sufficient for writing, in a straightforward way, a general-purpose 
NLO Monte Carlo algorithm. Note in particular that there is no need to extract a proper 
counter-term da"^ starting from a cumbersome expression for da^ in d dimensions. The 
NLO matrix element contributing to da^ can be evaluated directly in four space-time 
dimensions thus leading to an extreme simplification of the Lorentz algebra. 

The key point of our method for constructing the counter-term da^ is the dipole for- 
malism described in Sect. ^. Starting from the universal behaviour of the QCD matrix 
elements in the soft and collinear regions, we have introduced improved factorization for- 
mulae based on a dipole structure with respect to the momenta, colours and helicities of 
the QCD partons. Our dipole formulae correctly match the well-known singularities of 
the QCD scattering amplitudes in the soft and collinear limits. Moreover, these limits 
are approached smoothly, thus avoiding double counting of overlapping soft and collinear 
divergences. This smooth transition is possible because our dipole formalism is explicitly 
Lorentz covariant and the dipole formulae exactly fulfil momentum conservation. 

In the present paper, these main features of the dipole formulae have been used in 
the context of NLO computations of jet cross sections with no initial-state hadron. For 
lepton-hadron and hadron-hadron collisions, perturbative QCD calculations face additional 
difficulties related to the factorization of initial-state collinear singularities. Likewise, to 
calculate cross sections for processes in which a final state hadron is identified, the equiv- 
alent final-state collinear singularities must be dealt with. In a companion paper 
we show that the dipole formalism overcomes these difficulties in a simple and general 
manner. Thus we can explicitly provide a set of universal counter-terms or, equivalently, 
effective (non-singular) matrix elements that can be straightforwardly used for any NLO 
QCD calculation. 

^Actually, if the total number of QCD partons involved in the LO matrix element is less than or equal 
to three, one simply needs its incoherent sum over the colours (see Sect. ||). 
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